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We report on the status of a program of generating and using configurations with four flavors of 
dynamical quarks, using the HISQ action. We study the lattice spacing dependence of physical 
quantities in these simulations, using runs at several lattice spacings, but with the light quark mass 
held fixed at two tenths of the strange quark mass. We find that the lattice artifacts in the HISQ 
simulations are much smaller than those in the asqtad simulations at the same lattice spacings 
and quark masses. We also discuss methods for setting the scale, or assigning a lattice spacing to 
^ ensembles run at unphysical parameters. 



The XXVIII International Symposium on Lattice Field Theory, Lattice2010 
June 14-19, 2010 
Villasimius, Italy 



*Present address: Department of Physics, Brookhaven National Laboratory, Upton, NY 1 1973, USA 
' Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos. sissa. it/ 



Simulations with dynamical HISQ quarks 



D. Toussaint 3 



[000 



0.04 m„ 



-SB — 



OCCO 



CCCO 

0.045 



TOTO 



m 
o 
O 



txi 0.1 m s 

ccco 



[000 



TOTO 



O: generate 



O: potential O: f pi O: spectrum 



0.2 m s 
0.18 



TOTO 



TOTO 



0.15 



0.12 



4000 



TOTO 1 



0.09 



0.06 



HISQ program, April 6, 2010 

Figure 1: Status of the HISQ configuration generation program as of June 2010. We expect to eventually 
bring all these ensembles to approximately 1000 equilibrated configurations. The ensembles at the physical 
quark mass (0.04m 4 ) at a = 0.06 and 0.045 fm will require next generation computers. The black bars 
show the number of equilibrated configurations, and the red, blue and green bars show the progress of static 
potential, light pseudoscalar amplitudes, and simple hadron spectrum calculations. These are being done in 
parallel with the configuration generation to the extent possible. 



1. Introduction 



During the past eleven years the MILC collaboration has carried out a program of QCD sim- 
ulations using an improved staggered quark action, the "asqtad" action. This action has lattice 
artifacts that are much smaller than those in the original (one-link) staggered quark action. These 
simulations used two light and one strange dynamical quarks. Lattices have been generated with 
spacings ranging from 0.18 fm to 0.045 fm and light quark masses ranging from 0.05 m s to m s . The 
spatial sizes of these lattices ranged from 2.5 fm to 5.8 fm. In total, around 25000 configurations 
have been archived. These configurations have been used for a wide variety of QCD studies and 
are all publicly available from the NERSC Gauge Connection archive, from ILDG, or informally. 
A detailed review of this simulation program can be found in Ref. [1], 

Recently the HPQCD/UKQCD collaboration introduced a "highly improved staggered quark" 
(HISQ) action which further reduces taste symmetry violations, which are expected to be the largest 
lattice artifacts with staggered quark actions. Using this action, and taking advantage of improved 
machine power and algorithm development, we have begun a new program of lattice QCD simu- 
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lations. There are several differences between this new set of simulations and those in the asqtad 
program. First, the smeared links used in the D operator are done with a two level smearing, so 
that the quarks effectively see a smoother lattice, leading to reduced taste symmetry violations. 
Second, the three-link (Naik) term in D for the charm quark is modified to improve the disper- 
sion relation of the free charm quark. These two improvements together make up what is usually 
called the HISQ action[2]. In addition, this generation of simulations is using spatial lattice sizes 
about 20% larger than the comparable asqtad action simulations. While this seems like a small 
adjustment, since finite size effects are expected to decrease as exp(— m n V), these effects should be 
significantly reduced. Also, the Symanzik improved gauge action contains the one-loop effects of 
the HISQ fermions[3]. The comparable asqtad effects were not available when that program was 
started, but turned out to be surprisingly large[4]. With experience from the asqtad simulations, 
we expect to do a better job of tuning the dynamical quark masses, reaching a level of perhaps 2% 
accuracy. In contrast, some of the early asqtad lattice ensembles had strange quark masses as much 
as 20% off their a posteriori correct value. Finally, we are including a dynamical charm quark. At 
our smaller lattice spacings, am c is not large. We do expect the effects of dynamical charm to be 
small in most cases, but they may be important to the high temperature QCD part of our program, 
and they are quite cheap in include. 

In this program we expect to produce full scale ensembles with about 1000 equilibrated con- 
figurations at lattice spacings of 0.15, 0.12, 0.09 and 0.06 fm, with light quark masses of m s /5, 
m s /10 and the physical value, roughly m s /21. In the not too distant future we also expect to gener- 
ate an ensemble with a = 0.045 fm and the physical quark masses. We will also generate a number 
of coarser ensembles, mainly for setting the scale in HISQ high temperature QCD simulations. For 
each of these approximate lattice spacings we are using values of 10/g 2 and uo determined from the 
mi = m s /5 runs at all light quark masses. The strange and charm quark masses are adjusted from 
short tuning runs to give the physical values of 2M\ — M\ and \M\y + \M r[c respectively. Figure 1 
shows the current (June 2010) status of these runs. 

2. Tests of scaling 

The first stage in this program was to generate complete ensembles at a fixed but unphysical 
light quark mass of m s /5 at lattice spacings of 0.15, 0.12 and 0.09 fm. This allows us to test scal- 
ing, or dependence of calculated quantities on the lattice spacing. In particular, we wish to see if 
the reductions in taste symmetry violations are accompanied by reductions in lattice spacing de- 
pendence of other quantities. The results of these tests are briefly summarized here, but a complete 
report may be found in Ref. [5]. 

The HPQCD/UKQCD collaboration demonstrated the reduction of taste violations with HISQ 
valence quarks using quenched and asqtad sea quarks [2]. As expected, we find similar results with 
HISQ sea quarks. Figure 2 shows results for the taste splittings with HISQ valence and sea quarks 
with mi = ra. v /5, including some results on a partially completed ensemble with a = 0.06 fm. 

Perhaps the most stringent test of improvement is the topological susceptibility, since the de- 
pendence of this quantity on the sea quark mass demonstrates that the sea quarks are modifying 
the gluon configurations. Figure 3 shows the topological susceptibility for most of the asqtad 
ensembles[7], together with the HISQ results for mi = m s /5 at a = 0.12 and 0.09 fm. It can be seen 
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Figure 2: Taste splittings among the pions, comparing the asqtad and HISQ actions. The asqtad results used 
configurations with 2+1 flavors of dynamical quarks, and the HISQ results 2+1+1 flavors. The quantity plot- 
ted is r\ (m\ — Mq\ , where M n is the mass of the non-Goldstone pion and Mq is the mass of the Goldstone 
pion. Reading from top to bottom, the non-Goldstone pions are the % s (box), 7Co (fancy box), 7Tj (fancy plus), 
%io (plus), %ij (diamond), (cross) and 7To5 (octagon). r\ \M\—Mq\ is known to be almost independent 
of the light-quark mass. The vertical bar at the upper left shows the size of a factor of three, roughly the ob- 
served reduction in taste splittings, while the sloping solid line shows the theoretically expected dependence 
on lattice spacing. Nearly degenerate points have been shifted horizontally to improve their visibility. 



that the HISQ points fall below the corresponding asqtad points, which are indicated by the arrows 
in the figure. Note that the HISQ points are to the left of the corresponding asqtad points. This is 
because the horizontal axis is the mass of the taste singlet pion (the heaviest pion taste), and the 
reduction in taste symmetry breaking moves the points to the left. It is the movement down relative 
to the asqtad points that represents an improvement in the gluon configurations. 

We also looked at the masses of the p and nucleon in these ensembles. Strictly, these are an 
unphysical vector meson and baryon with both valence and sea light quark masses at m s /5. Thus, 
we can check their dependence on lattice spacing, but without a chiral extrapolation we cannot test 
their masses against experiment. With this caution, we found that the dependence of r\M p and 
r\Mf] on the lattice spacing was much smaller than for the asqtad ensembles at comparable quark 
masses. More details are in Ref. [5]. Of course, this could be interpreted either as improved scaling 
of hadron masses using r\ to set the scale or improved scaling of n using hadron masses to set the 
scale. 
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Figure 3: The topological susceptibility, comparing the asqtad and HISQ actions. Points with the asqtad 
action are shown for several lattice spacings and quark masses, and the HISQ results for a ps 0.12 fm and 
a ps 0.09 fm with m\ = 0.2m,. For the horizontal axis we use the mass of the taste singlet pion, since in lowest 
order chiral perturbation theory the topological susceptibility is a function of this mass [6]. The curves in 
the figure come from a chiral perturbation theory fit to the asqtad data. A discussion of the methods used 
for computing the susceptibility and the asqtad results can be found in Ref. [7]. The two arrows indicate the 
locations of asqtad points with lattice spacing and quark mass similar to the two HISQ points. (In the case 
of the a ps 0.09 fm HISQ point, the quark mass falls between two of the masses of the asqtad points.) 

3. Setting the length scale 

QCD simulation programs involving unphysical quark masses generally involve a definition 
of the lattice spacing as a function of the gauge couplings and sea quark masses. About the only 
real requirement on this definition is that it should be correct at the physical quark masses in the 
continuum limit. The most common choice is a scale determined from the static quark potential[8], 
usually either r\F{jo) = —1.65 or r\F(r\) = — 1. This is not a physical quantity, but it is a conve- 
nient interpolating quantity that can be determined accurately with a reasonable amount of work. 
Its value is determined by matching some physical quantity, such as quarkonium mass splittings 
or f % , in the continuum and chiral limits. Recently the HPQCD collaboration has suggested using 
the decay constant of a fictitious isovector pseudoscalar meson with valence quark masses equal to 
the strange quark mass as the length standard[9], which we will call f ss . Like ro or r\, this is an 
unphysical quantity whose value is determined by matching a physical quantity, most likely f n , in 
the continuum limit at the physical quark mass. Of course, one could use any mass for the valence 
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Figure 4: Differences in determinations of the length scale using different standards. In the legend, the 
symbol types are labelled as "valence on sea". The "FP4S" points (bursts) use a pseudoscalar meson with 
both valence quarks at 0.4 times the strange quark mass, while all the other points use valence quarks at the 
strange quark mass. The "HISQ on asqtad" points are taken from Ref. [9]. 



quark, and we have been experimenting with using a mass of 0.4 times the strange quark mass. The 
choice of reference mass involves a tradeoff — we want a mass heavy enough so that an accurate 
simulation result is possible but light enough so that chiral perturbation theory can be used with 
confidence. 

Use of an unphysical decay constant has several advantages. It is possible to get very good 
accuracy, and the systematic errors from things like choices of fit ranges are better understood than 
those in r\. Also, for calculations involving hadron masses and matrix elements, there are fewer 
steps in the logic relating the lattice results to dimensionful quantities. On the other hand, the meson 
decay constant takes much more computer time than the static quark potential (although it may be 
needed anyway) and is dependent on the choice of valence quark formulations. For example, 
asqtad and HISQ valence quarks would give different lattice spacings for the same ensemble of 
configurations. 

Table 1 shows the lattice spacings for several asqtad and HISQ ensembles determined from 
H and from pseudoscalar amplitudes with asqtad and HISQ valence quarks. Figure 4 shows the 
differences between the pseudoscalar and r\ lattice spacings versus a 2 , including lattice spacings 
determined from both /„ and f ^ s , and it can be seen that the differences are vanishing in the 
continuum limit. 
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Action 10/g 2 am/ am s am c 


a{r\) a(/M,asq) o(/m,hisq) 


asqtad 6.76 0.01 0.05 
asqtad 7.09 0.0062 0.031 
asqtad 7.46 0.0036 0.018 
asqtad 7.81 0.0028 0.014 


0.1178(2) 0.1373(2) 0.1264(11) 
0.0845(1) 0.0905(3) 0.0878(7) 
0.0588(2) 0.0607(1) 0.0601(5) 
0.0436(2) 0.0444(1) 0.0443(4) 


HISQ 5.80 0.013 0.065 0.838 
HISQ 6.00 0.0102 0.0509 0.635 
HISQ 6.30 0.0074 0.037 0.440 


0.1527(7) na 0.1558(3) 
0.1211(2) na 0.1244(2) 
0.0884(2) na 0.0900(1) 



Table 1: Lattice spacing determinations by various methods. The values for HISQ valence quarks with 
asqtad sea quarks are taken from Ref. [9]. In the other cases, the errors are statistical only; they do not 
include the errors on the physical value of r\ (n = 0.3108(15)(+^)) or on the estimate of f ss in MeV, 
/„ = 181.5(1.0) MeV[9]. 
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